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Current phenomenological studies of jet observables at colliders are clearly limited by the theoretical uncer- 
tainties inherent in the next-to-leading order QCD description. We discuss the recent progress made towards 
the calculation of QCD corrections to jet observables at the next-to-next-to-leading order in QCD and highlight 
future perspectives and yet open issues. 



1. Introduction 

Jet production observables are among the most 
sensitive probes of QCD at high energy colliders, 
where they are used for example to determine the 
strong coupling constant. At present, the inter- 
pretation of jet production data within pertur- 
bative QCD is restricted to next-to-leading order 
(NLO) calculations, with theoretical uncertain- 
ties considerably larger than current experimen- 
tal errors. Going beyond NLO calculations offers 
moreover a more accurate matching of theoreti- 
cal and experimental jet definitions and a more 
detailed modelling of the hadronic final state [ 
[[). The extension of jet calculations to NNLO 
requires three ingredients: the two-loop correc- 
tions to multi-leg amplitudes, the single unre- 
solved limits of one-loop amplitudes and the dou- 
ble unresolved limits of tree amplitudes. Finally, 
all these contributions have to be combined to- 
gether into a program for the numerical com- 
putation of jet observables from NNLO parton 
level cross sections. In this talk, I review recent 
progress made on these subjects as well as the 
currently open issues. 

2. Virtual Two-Loop Corrections 

Within dimensional regularization, the large 
number of different integrals appearing in multi- 
loop calculations can be reduced to a small 
number of so-called master integrals by us- 
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ing integration-by-parts (IBP) identities [ 0, [|. 
These identities exploit the fact that the integral 
over the total derivative of any of the loop mo- 
menta vanishes in dimensional regularization. 

For integrals involving more than two exter- 
nal legs, another class of identities exists due 
to Lorentz invariance. These Lorentz invariance 
identities [ Q rely on the fact that an infinites- 
imal Lorentz transformation commutes with the 
loop integrations, thus relating different integrals. 
The common origin of IBP and LI identities is 
the Poincare invariance of loop integrals within 
dimensional regularization, as was pointed out 
by J.J. van der Bij at this conference. Using 
integration-by-parts and Lorentz invariance iden- 
tities, all two-loop Feynman amplitudes for 2 — > 2 
scattering or I — > 3 decay processes can be ex- 
pressed as linear combinations of a small number 
of master integrals, which have to be computed 
by some different method. Explicit reduction for- 
mulae for on-shell two-loop four-point integrals 
were derived in [ |j| . Computer algorithms for the 
automatic reduction of all two-loop four-point in- 
tegrals were described in [ |[ |J. 

The master integrals relevant to 2 — > 2 scatter- 
ing or 1 — > 3 decay processes are massless, scalar 
two-loop four-point functions with all legs on- 
shell or a single leg off-shell. Several techniques 
for the computation of those functions have been 
proposed in the literature, such as the application 
of a Mellin-Barnes transformation to all propaga- 
tors, which was used successfully to compute the 
on-shell planar double box integral [ ^, 
on-shell non-planar double box integral [ 
two double box integrals with one leg off-shell 
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Figure 1. Computer Algebra for analytic evaluation of two- loop matrix elements. 



[T(j]. Most recently, the same method was used 
to derive the on-shell planar double box integral 
with one internal mass scale [ as well as the 
high energy limit of the on-shell planar triple box 
integral [ 

A method for the analytic computation of mas- 
ter integrals avoiding the explicit integration over 
the loop momenta is to derive differential equa- 
tions in internal propagator masses or in exter- 
nal momenta for the master integral, and to solve 
these with appropriate boundary conditions. The 
computation of master integrals from differential 
equations proceeds as follows [ |ij . Carrying out 
the derivative with respect to an external invari- 
ant on the master integral of a given topology, 
one obtains a linear combination of a number of 
more complicated integrals, which can however 
be reduced to the master integral itself plus sim- 
pler integrals by applying the reduction methods 



discussed above. As a result, one obtains an inho- 
mogeneous linear first order differential equation 
in each invariant for the master integral. The in- 
homogeneous term in these differential equations 
contains only topologies simpler than the topol- 
ogy under consideration, which are considered to 
be known if working in a bottom-up approach. 
The master integral is then obtained by match- 
ing the general solution of its differential equation 
to an appropriate boundary condition. 

Using the differential equation technique, one 
of the on-shell planar double box integrals [ 
as well as the full set of planar and non-planar 
off-shell double box integrals [ were derived. 
The computer algebra structures applied in the 
computation of the master integrals from the dif- 
ferential equations are displayed in the right hand 
column of Figure |l|. The differential equation ap- 
proach has recently been extended to phase space 
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integrals [ |l5| . 

A strong check on all these computations of 
master integrals is given by the completely nu- 
merical calculations of [ [l6| , which are based on 
an iterated sector decomposition to isolate the in- 
frared pole structure. The methods of [ [ll| were 
applied to confirm all of the above-mentioned cal- 
culations. 

A third approach, which avoids the reduction 
to master integrals, has been presented in [ p7[ . 
In this approach, all integrals appearing in the 
two-loop amplitudes are related to higher tran- 
scendental functions, which can be expanded in 
terms of nested harmonic sums. 

The two-loop four-point functions with all legs 
on-shell can be expressed in terms of Nielsen's 
polylogarithms [ ^9|. In contrast, the closed 
analytic expressions for two-loop four-point func- 
tions with one leg off-shell contain two new classes 
of functions: harmonic polylogarithms [ 
and two-dimensional harmonic polylogarithms 
(2dHPL's) [ Accurate numerical implemen- 
tations for these functions [ are available. 

2.1. 2 — > 2 Processes with all legs on-shell 

With the explicit solutions of the integration- 
by-parts and Lorentz-invariancc identities for on- 
shell two-loop four-point functions and the corre- 
sponding master integrals, all necessary ingredi- 
ents for the computation of two-loop corrections 
to 2 — > 2 processes with all legs on-shell are now 
available. The generic structure of such a cal- 
culation is outlined in the left hand column of 
Figure H. In fact, only half a year elapsed be- 
tween the completion of the full set of master in- 
tegrals and the calculation of the two-loop QED 
corrections to Bhabha-scattering [ ^2). Subse- 
quently, results were obtained for the two-loop 
QCD corrections to all parton-parton scattering 
processes [ p3| . For gluon-gluon scattering, the 
two-loop helicity amplitudes have also been de- 
rived [ Moreover, two- loop corrections were 
derived to processes involving two partons and 
two real photons [ [25). Since the gluon fusion 
into photons has a vanishing tree level amplitude, 
these results form part of the NLO corrections to 
photon pair production [ [2(|, yielding a sizable 
correction. 



Finally, light-by-light scattering in two-loop 
QED and QCD was considered in [ H?J . 

The results for the two-loop QED matrix ele- 
ment for Bhabha scattering [ ^] were used in [ 
p8[ to extract the single logarithmic contributions 
to the Bhabha scattering cross section. 

2.2. 2^2 Processes with one off-shell leg 

Using the two-loop master integrals with one 
off-shell leg [ [14|, the two-loop QCD matrix ele- 



ment for 



3 jets 



29[ and the correspond- 
ing helicity amplitudes [ tMI were computed fol- 
lowing the reduction procedure depicted in Fig- 
ure |l|. The infrared pole structure of these results 
agrees with the prediction [ [}lj obtained from an 
infrared factorization formula. 

An independent confirmation of part of these 
results was performed in [ |3^] , where two of the 
seven colour factors (corresponding to the terms 
proportional to nj) of the two- loop helicity am- 



plitudes for e" 1 



3 jets were derived using the 



nested sum method of [ |l7|] . 

Processes related to e + e~ — > 3 jets by cross- 
ing symmetry are (2 + l)-jet production in deep 
inelastic ep scattering and vector-boson-plus-jet 
production at hadron colliders. The analytic con- 
tinuation of the e + e~ — > 3 jets two-loop helicity 
amplitudes to the kinematic regions relevant for 
these scattering processes has been derived in [ 

E- 

3. Real Corrections 

Besides the two- loop virtual corrections, a 
NNLO calculation of jet observables has to in- 
clude the contributions from single unresolved 
(soft or collinear) real radiation from one-loop 
processes as well as from double unresolved real 
radiation at tree level. Only after summing all 
these contributions (and including terms from the 
renormalization of parton distributions for pro- 
cesses with partons in the initial state), do the 
divergent terms cancel among one another. The 
factorization properties of both the one-loop, one- 
unresolved-parton contribution [ B4| and the tree- 



level, two-unresolved-parton contributions [ 35 
have been studied, but a systematic procedure for 
isolating the infrared singularities has so far been 
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Figure 2. Structure of NNLO parton level Monte Carlo programme. 



established only for the one-loop, one-unresolved- 
parton processes. Although this is still an open 
and highly non-trivial issue, significant progress 
is anticipated in the near future. 

4. Numerical Implementation 

Prior to their implementation into a numer- 
ical program, it is needed to analytically ex- 
tract the infrared pole terms from the one- 
loop, one-unresolved-parton and two-unresolved- 
parton contributions. At NLO, two types of 
methods have been used very successfully in the 
past. The phase space slicing method [ [36) divides 
up the final state phase space into resolved and 
unresolved regions; the infrared subtraction terms 
are then integrated only over the unresolved re- 
gions. This method avoids overcounting of singu- 



lar contributions, the required integrals over re- 
stricted phase space regions can however be very 
involved beyond NLO. The subtraction method [ 
|37j integrates the subtraction terms over the full 
phase space. The construction of the subtraction 
term requires in this case great care to avoid over- 
counting problems. In general, algorithms based 
on subtraction are more efficient numerically. 

The remaining finite terms must then be com- 
bined into a numerical program implementing 
the experimental definition of jet observables and 
event-shape variables. The sketch of such a pro- 
gramme to compute e + e~ — > 3j at NNLO is given 
in Figure [|. 

Programs to compute processes with initial 
state hadrons involve the additional complica- 
tion of initial state singularities, which have to 
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be absorbed into the NNLO parton distributions. 
Lacking the full expressions for the splitting func- 
tions at this order, these are not yet available at 
present, work on them is however well advanced [ 



A first calculation involving the features of 
NNLO jet calculations was presented for the case 
of photon-plus-one-jet final states in electron- 
positron annihilation in [ , thus demonstrating 
the feasibility of this type of calculations. A pre- 
requisite for such a numerical program computing 
n jet final states is a stable and efficient next-to- 
leading order programme for the processes yield- 
ing n + 1 jet final states. For the processes of 
highest phenomenological interest, these are al- 
ready available: e + e~ — ► 4j [ [jo], ep — » (3 + l)j [ 

5. Conclusions and Outlook 

Considerable progress has been made in the 
last two years (in fact since the last "Loops and 
Legs" - and "RADCOR" -conferences) towards the 
computation of jet observables at NNLO in QCD. 
In particular, new methods have been developed 
for the calculation of two-loop virtual corrections 
to four-point scattering amplitudes. As a result, 
the two-loop virtual corrections relevant to all 
phcnomcnologically important processes in QCD 
and QED are now known. 

These form however only part of the full calcu- 
lation required at NNLO accuracy, which also has 
to take into account contributions from one-loop 
single-unresolved radiation and tree-level double 
unresolved radiation processes. While appropri- 
ate subtraction terms for these processes have 
been known for quite some time, their analytic 
integration (required for the cancellation of in- 
frared poles in physical jet observables) is still an 
unsolved problem. Once this obstacle has been 
overcome, the remaining finite parts can be im- 
plemented into a numerical programme to com- 
pute NNLO jet observables. 
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